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Abstract 

^ Nernst effect, the transverse voltage generated by a longitudinal thermal gradient 

in presence of magnetic field hjis recently emerged as a very sensitive, yet poorly 

^ understood, probe of electron organization in solids. Here we report on an experiment 

^ on graphite, a macroscopic stack of graphene layers, which establishes a fundamental 

^ link between dimensionality of an electronic system and its Nernst response. In 

I 

sharp contrsist with single-layer graphene, the Nernst signal sharply peaks whenever 

2 a Landau level meets the Fermi level. This points to the degrees of freedom provided 

by finite interlayer coupling as a source of enhanced thermoelectric response in the 

vicinity of the quantum limit. Since Landau quantization slices a three-dimensional 

^ Fermi surface, each intersection of a Landau level with the Fermi level modifies the 

^ Fermi surface topology. According to our results, the most prominent signature of 

On 

O such a topological phase transition emerges in the transverse thermoelectric response. 

G> 

O 
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Graphene, a single layer of carbon atoms set in a honeycomb lattice, has attracted much 
attention because of its uncommon properties [T]. The two-dimensional gas of massless elec- 
trons embedded in graphene displays a particular version of Quantum Hall Effect [21 E]- Its 
thermoelectric response has been a subject of recent experimental investigations by several 
groups in El E], who found results in agreement with the theory of thermoelectricity for a 
two-dimensional electron system in the quantum Hall regime[71 [8]. On the other hand, the 
thermoelectric tensor of graphite [9j, a macroscopic stack of graphene layers, has never been 
studied at low temperatures and high magnetic fields. Here we report on such a study un- 
covering the crucial role of interlayer coupling. When the field is strong enough to push the 
system to the quantum limit, the quantum oscillations of the Nernst response easily dwarf 
the oscillations seen in other transport coefficients including the Seebeck signal. Moreover, 
the Nernst signal sharply peaks whenever a Landau level intersects the chemical potential. 
Both these features are absent in graphene, but were previously reported in bulk bismuth 
near the quantum limit |1U]. The resolved Nernst peaks emerge as signatures of a topological 
phase transition [TT], which occurs at the intersection of Landau and Fermi level [12]. Such a 
field-induced modification of the Fermi surface topology at lower Landau levels is exclusive 
to three dimensions, for which no adequate quantitative description of the thermoelectric 
response in the vicinity of the quantum limit is yet available. The results point to the con- 
figurational degrees of freedom associated with the finite dispersion of electrons along the 
magnetic field as the source of a huge off-diagonal thermoelectric response. 

Fig. 1 presents the thermal evolution of the field dependence of the Nernst signal, S^;^ 
in two highly-oriented pyrolytic graphite(HOPG) samples. As the temperature decreases, 
quantum oscillations become sharper and their visibility extends to lower fields. In graphite, 
quantum oscillations of various physical properties such as resistivity (the Shubnikov-de Hass 
effect )[T3]. magnetic susceptibility (the de Hass- van Alphen effect) [H], the Hall coefficient 
and thermopower[T7] were measured many years ago and the results are in agreement with 
the structure of the Fermi surface predicted by the Slonczewski-Weiss-McClure (SWM) 
model[T5l [16] . More recently, quantum oscillation studies employing new techniques for 
analysis and measurements have emerged and the degree of accuracy by which the SWM 
model describes the fine structure of the Fermi surface has become a subject of ongoing 
debate [TSl[TU] . However, in the case of the Nernst effect, the only available report is restricted 
to temperatures above the liquid nitrogen [2U]. 
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The particularity of the Nernst effect as a probe of quantum oscillations at low index Lan- 
dau levels is highlighted in Fig. 2. Quantum oscillations are visible in the field dependence 
of both diagonal and off-diagonal components of electric and thermoelectric conductivity 
tensors. However, only in the case of transverse thermoelectric response, the oscillating 
signal peaks by far the non-oscillating background. In contrast, the Shubnikov-de Haas 
oscillations represent a tiny fraction of the overall longitudinal resistivity. Moreover, while 
the longitudinal response dominates for resistivity {p^x ^ Pxy)i the Nernst signal is an order 
of magnitude larger than the Seebeck coefficient {S^y ^ S^x)- The large magnitude of the 
Nernst coefficient in graphite (second only to bismuth|21j among metals) is an expected 
consequence of the high mobility (3 x l^f' crn^V^^ s^^) and the low Fermi energy (19 meV) 
of the system (See the supplementary information for details). 

Fig. 3 compares the fine structure of pxy and S^y when the lowest Landau levels cross 
the Fermi level. The Nernst response peaks when the magnetic field attains a value corre- 
sponding to an intersection of a Landau level with the the Fermi level. The magnitude of 
these fields were determined by theory [22] and confirmed by experiment [T7] many years ago 
(See supplementary information for details on the indexing of the Landau levels and a com- 
parison of the fine structure in Sxy and S^^;)- The peaks in S^^j, for the lowest Landau levels 
of the two pockets can be clearly associated with a jump or drop in p^y depending on the 
sign (hole-like or electron-like) of the carriers of the pocket. Between jumps and drops, that 
is when the Fermi level lies between distinct Landau levels, the Hall resistivity, while not 
strictly fiat, shows little field dependence. Presenting the same data in a different fashion 
in the lower panel of the same figure, one can see that at high fields and low temperatures, 
|; becomes constant between two successive peaks. In other words, in the low-T high-B 
limit, when a Landau level does not intersect the Fermi level , the Nernst response becomes 
both 5-linear and T-linear but sharply peaks otherwise, with a singularity which becomes 
more pronounced with cooling. In addition to the two above-mentioned HOPG samples, we 
also measured a natural single crystal of graphite and found very similar results (see the 
supplementary information). 

It is illuminating to compare these features with results very recently reported for 
graphene in E] . The first difference is the sheer magnitude of the Nernst response. Whereas 
in graphene, the measurable signal at T=10 K is in the range of 10 to 20 pVK~^, the 
Nernst signal in the same temperature range in graphite is two orders of magnitude larger 
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and approaches 1 mV (Fig. la). This difference is most probably due to a difference 
in electron mobility. In graphite, it can exceed 10^cm^y~^s~^[23], two orders of magnitude 
larger than most currently available graphene samples, including those used in the Nernst 
studies. The second and more fundamental difference is the field profile of the Nernst re- 
sponse. In graphene, the intersection of the chemical potential and a Landau level (which 
can be realized either by scanning the field or modifying the gate voltage) leads to a change 
in the sign of S^.^^. In other words, the coincidence of a Landau level with the chemical poten- 
tial is concomitant with a jump in p^y and a vanishing S^y In graphite, as seen above, the 
Nernst response attains its maximum in the same conditions. Note that the huge electron 
mobility, the origin of the large low-field non-oscillating Nernst response in graphite, can 
neither explain the preponderance of the oscillating component, nor its particular profile. 

Interestingly, the field profile of the Nernst response we find here is analogous to the 
one previously reported in bulk bismuth |10j. On the other hand, the Nernst response of 
graphene in E] presents a functional form reminiscent of the case of the two-dimensional 
electron systems realized in semiconductor heterostructures [25] . In both latter cases, '^^y 
vanishes at the intersection of Landau and Fermi levels as predicted by the theory conceived 
for two dimensions O El l2l] . 

Two features set apart graphene from both bismuth and bulk graphite. The most fun- 
damental is the presence of a finite third-axis dispersion in the bulk materials. A finite 
interlayer coupling, no matter how small, would warp the perfect cylindrical Fermi surface 
of a two-dimensional monolayer. According to the SWM model, the Fermi surface of both 
electrons and holes in graphite are elongated ellipsoids. The precise magnitude of inter- 
layer coupling has been a matter of debate and the Fermi surface of one type of carriers 
may not be closed along the c-axis. In any case, however, it cannot be a perfect cylinder. 
The second difference is that both bismuth and bulk graphite are compensated semi-metals 
hosting equal concentrations of mobile carriers of both sign. Experiments on graphene were 
performed on a system with a single type of carriers determined by the sign of the applied 
gate voltage. These two features lead to the emergence of a qualitatively different transverse 
thermoelectric response. 

The consequences of the first difference is sketched in fig. 4, which compares the passage 
of successive Landau levels through the chemical potential in presence (3D, left ) or absence 
(2D, right) of the z-axis dispersion. In three dimensions. Landau quantization truncates the 
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Fermi surface. The wave-vector of mobile electrons are confined to slices formed by the inter- 
section of the Fermi-Dirac distribution and Landau spectrum. The thickness of these slices 
depends on temperature and disorder. As the field is swept the slices move outward and 
downward. The Nernst peaks are concomitant with the merger of two slices. As first noticed 
many years ago[12j. such a merger is a case of electronic topological phase transition [Hj. In- 
terestingly, a singularity in the Nernst response in the vicinity of a topological transition was 
theoretically predicted in another context [26]. On the other hand, in the two-dimensional 
case, sketched in the right side of the same figure, no such topological transition occurs when 
the ellipsoid is replaced by a perfect cylinder. Here, when the Landau level coincides with 
chemical potential, the system becomes dissipative, otherwise it is gapped. 

In the two-dimensional case, displacing a Landau level across the chemical potential 
leads to a symmetric Sxy profile. This can be qualitatively understood in the semi-classical 
picture. The thermoelectric tensor S measures the change in conductivity tensor, a, caused 
by a small shift in the chemical potential: 

~ TT^klT da 

According to this equation (sometimes called the extended Mott formula), if the chemical 
potential happens to be at its optimal position for a maximal Hall mobility the transverse 
thermoelectric response, S^j^, is expected to vanish. This simple argument can explain the 
qualitative profile of the Nernst response in the 2D case, which is symmetric and vanishes at 
the intersection. Since when a Landau level is precisely at the chemical potential, any shift 
of the latter would diminish the Hall mobility. This zero response should be sandwiched 
between a positive and a negative peak (as sketched in Fig.4d) . The sign of these two reflect 
the sign of the energy derivative when one approaches the critical field from lower fields or 
higher fields. 

The profile of the Seebeck coefficient in the 3D graphite (See Fig. S3 in the supplementary 
information) is not qualitatively different from the 2D one. On the other hand, the Nernst 
response in the 3D configuration emerges as a remarkable singularity. The clear asymmetry 
between configuration a and configuration c in Fig. 4 is the key to the functional shape 
of the Nernst response. Experimentally, Sxy steadily increases as the two separated slices 
of the Fermi surface move towards each other. It peaks when they merge to form a single 
one and begins a sudden drop after the Landau level moves beyond the Fermi level. The 



configurational entropy of the electronic states with vanishing and quantized k^y, which 
emerge as this topological phase transition occurs, appears to be the source of the enhanced 
Nernst response exclusive to the 3D case. 

The compensated nature of the bulk semi-metals under study imposes another constraint. 
Charge neutrality implies that any change in the density of carriers with one sign would 
shift the chemical potential in a way to maintain the equality between the concentration of 
electron and holes. A topological transition in one (say hole-like) ellipsoid thus moves the 
chemical potential and modifies the shape of the other. 

When electrons are confined to two dimensions, the quantum oscillations of the Hall 
response become sharp steps in the vicinity of the quantum limit. It emerges from this study 
that such a confinement modifies also the quantum oscillations of the Nernst response in a 
remarkable way.A satisfactory theoretical understanding of the Nernst effect in bulk systems 
across the quantum limit is still missing. It is a remarkable irony that the least theoretically 
understood transport coefficient happens to be the most experimentally sensitive probe of 
quantum oscillations. 

Extending theses measurements to higher magnetic fields would let one probe the ther- 
moelectric response of graphite in the extreme quantum limit. Such a study in bismuth 
has recently uncovered several enigmatic field scales[27]. One collateral conclusion of study 
presented here is to rule out that the Nernst anomalies seen there could originate from the 
two-dimensional surface states. Indeed, the field profile of the Nernst anomalies observed in 
bismuth clearly points to a three-dimensional origin. 
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FIG. 1: Zhu et al. 



Fig. 1. The field dependence of tlie Nernst signal in two Highly-Oriented Pyrolytic Graphite 
(HOPG) samples at different temperatures. The signal was measured with an in-plane 
thermal gradient and a field oriented perpendicular to the layers. As the temperature is 
lowered, oscillations become sharper and the Zeeman splitting becomes visible. 



9 



3 6 9 12 




3 6 9 12 



B(T) 

FIG. 2: Zhu et al. 

Fig. 2. Quantum oscillations of various transport coefficients at low temperatures. The up- 
per panel shows longitudinal and Hall resistivity of sample 2. Shubnikov-de Haas oscillations 
superpose on a huge non-oscillating background. The lower panel shows the field depen- 
dence of the thermoelectric coefficients in the same sample. The Nernst response dominates 
by far the Seebeck coefficient and presents pronounced well-defined oscillations. The Fermi 
surface and the Brillouin zone of graphite are shown as an inset in the upper panel. The 
Fermi surface consists of 6 pairs of adjacent ellipsoid pockets (in blue and green) hosting 
carriers of opposite signs. 
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FIG. 3: Zhu et al. 

Fig. 3. The variation witli inverse magnetic field of tlie Hall resistivity,p2.j^ (a), the Nernst 
signal, S^y (b), and the Nernst coefficient divided by temperature, ^ = at 3 different 
temperatures (c). When a Landau level meets the Fermi level, the Hall resistivity presents a 
jump and the Nernst signal sharply peaks. At low temperature and high field, Sxy becomes 
linear in field and temperature when the Fermi level is between two Landau levels. The 
inset zooms at peak positions. When a Landau level is at the Fermi level, ^ increases with 
decreasing temperature, otherwise it attains a temperature-independent magnitude. Landau 
levels of the pockets are indexed in panel b. 
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FIG. 4: Zhu et al. 

Fig. 4. A comparison of liow Landau levels (LLs, indexed n, n-1) cross the chemical 
potential, /x, as the magnetic field is swept, in the 3D (panels a to c) and 2D (panels a' to 
c') cases. In each panel the Fermi surface with cuts of the Landau cylinders in reciprocal 
space is shown together with a plot of the corresponding energies. The z-axis is set by 
the direction of the applied magnetic field. In the 3D case , an ellipsoidal Fermi surface 
is truncated to slices, which are the intersections of the Fermi level and the Landau levels. 
In the 2D case, the Fermi surface is a perfect cylinder and there are no slices. The Nernst 
signal is sketched in the panels d and d'. In 3D, peaks when the Landau level meets 
the Fermi level and a topological phase transition occurs by the merger of two slices. The 
odd number of slices in panel b is to be contrasted with the even number of slices in panel 
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a and c. In 2D, Sxy vanishes when a Landau level intersects the chemical potential. The 
experimentally resolved Nernst signal presents the profile of panel d in both graphite [as 
reported here] and bismuth |10|. On the other hand, in graphenejH E] and semiconducting 
heterostructures|25j it presents the profile of panel d'. 
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